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Abstract 

Quantum interference is investigated within the complex quantum Hamilton- Jacobi formalism. As shown in a previous work [Phys. 
H Lett. 102, 250401 (2009)], complex quantum trajectories display helical wrapping around stagnation tubes and hyperbolic 
deflection near vortical tubes, these structures being prominent features of quantum caves in space-time Argand plots. Here, we 
C^l further analyze the divergence and vorticity of the quantum momentum function along streamlines near poles, showing the intricacy 
^^of the complex dynamics. Nevertheless, despite this behavior, we show that the appearance of the well-known interference features 
(on the real axis) can be easily understood in terms of the rotation of the nodal line in the complex plane. This ofl'ers a unified 
^ description of interference as well as an elegant and practical method to compute the lifetime for interference features, defined in 
terms of the average wrapping time, i.e., considering such features as a resonant process. 

CN 'Keywords: Quantum interference. Complex quantum trajectory. Quantum momentum function, Polya vector field. Quantum cave 
'PACS: 03.65.Nk, 03.65.Ta, 03.65.Ca 
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^ X' Introduction 

a ■ 

d Quantum interference is an observable eff'ect arising from 
^^the coherent superposition of quantum probability amplitudes. 
, ^^ It is involved in a very wide range of important applications, 
jsuch as superconducting quantum interference devices or SQUIDs 
ill, 01], coherent control of chemical reactions |3], atomic and 
molecular interferometry |4-7], and Talbot/Talbot-Lau inter- 
ferometry with relatively heavy particles (e.g., Na atoms (sl) 
and Bose-Einstein condensates |9]). Indeed, possibly one of 
the main practical applications of interference nowadays is in 
Bose-Einstein condensate (BEC) interferometry Ill0l4l21 due to 
its potential use in applications, such as sensing, metrology or 
Q\ quantum information processing. Thus, since the first experi- 
mental evidence of interference between two freely expanding 
^ BECs was observed 1 13], an increasing amount of work, both 

theoretical and experimental, can be found in the literature 
r \ [2O]. Basically, in this type of interferometry there are three 
d steps. First, the atomic cloud is cooled in a magnetic trap until 
condensation takes place; second, the BEC is split coherently 
by means of radio-frequency 1 18] or microwave fields 1 19]; and 
third, the double- well-like trapping potential is switched off' and 
the two parts of the BEC are allowed to interfere by free ex- 
pansion (see, for example. Fig. 6 in Ref. |20] as an illustration 
of a real experimental outcome). Apart from the practical ap- 
plications mentioned before, interference plays also an impor- 
tant role when dealing with multipartite entangled systems as 
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an indicator of the loss of coherence induced by the interac- 
tion between the diff'erent subsystems, commonly referred as 
"Schrodinger cat states" | 21]. However, very little attention be- 
yond the implications of the superposition principle has been 
devoted to understanding quantum interference at a more fun- 
damental level. 

Bohmian mechanicsprovides an alternative interpretation 
to quantum mechanics [22142411 . As an analytical approach, this 
formulation has been used to analyze atom- surface scattering 
the quantum Talbot eff'ect 1I28I1 . quantum nonlocality 
1 29] or quantum interference |[3Q|i . As a synthetic approach, 
the quantum trajectory method (QTM) |[3lll has been developed 
as a computational implementation to the hydrodynamic for- 
mulation of quantum mechanics to generate the wave function 
by evolving ensembles of quantum trajectories | 32]. However, 
computational difficulties resulting from interference eff'ects are 
encountered in regions where the wave packet is reffected from 
the barrier and nodes or quasi-nodes occur. Thus, a bipolar 
counter-propagating decomposition approach for the total wave 
function has been developed to overcome numerical instabili- 
ties due to interferences 113314391] . 

Here, a more detailed analysis of quantum interference within 
the complex quantum Hamilton- Jacobi formalism [40, 41] is 
presented. The complex quantum trajectory method has been 
applied to stationary problems 1^-58] and to one-dimensional 
and multi-dimensional wave packet scattering problems 
[59I1 . The purpose of this work is to investigate how information 
can be extracted from this complex formulation regarding inter- 
ference. Recently, quantum interference associated with super- 
positions of Gaussian wave packets was thoroughly analyzed 
using both the standard quantum-mechanical and Bohmian ap- 
proaches in real space Ii3^i6fl]. In spite of its simplicity, such 
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superpositions are experimentally realizable in atom interfer- 
ometry |[6ll - [65ll , for example. In contrast to conventional quan- 
tum mechanics, Bohmian mechanics offers a trajectory-based 
understanding of quantum interference. The interference of 
two wave packets in one real coordinate leads to the forma- 
tion of nodal structure, and the quantum potential near these 
nodes forces these trajectories to avoid these regions and to ex- 
hibit laminar flow in space-time plots In contrast, within 
the complex quantum Hamilton- Jacobi formalism, the collinear 
collision of two Gaussian wave packets demonstrated caustics 
and vortical dynamics in the complex plane ||63,|6l]. This com- 
plicated trajectory dynamics cannot be found in Bohmian me- 
chanics unless two or more real coordinates are introduced. 

When systematically analyzed for local structures of the 
quantum momentum function (QMF) and the Polya vector field 
(PVF) around its characteristic points ||67l|680, complex quan- 
tum trajectories display helical wrapping around stagnation tubes 
and hyperbolic deflection around vortical tubes |66]. This in- 
triguing topological structure is formed by these tubes and gives 
rise to quantum caves in space-time Argand plots. Quantum in- 
terference thus leads to the formation of quantum caves and 
the appearance of the topological structure mentioned before. 
In contrast with quantum trajectories, Polya trajectories display 
hyperbolic deflection around stagnation tubes and helical wrap- 
ping around vortical tubes. The QMF divergence and vorticity 
characterize the turbulent flow of trajectories, determining the 
so-called wrapping time |66] for an individual trajectory. More- 
over, it is shown that both the PVF divergence and vorticity van- 
ish except at poles, thus the PVF describing an incompressible 
and irrotational flow. 

Trajectories launched from difl'erent positions wrap around 
the same stagnation tubes. Hence, the circulation of trajectories 
can be viewed as a resonance process, from which one can ob- 
tain a natural way to define a "lifetime" for the interference. 
This information could be therefore used with practical pur- 
poses to analyze, explain and understand experiments where 
interference is the main physical process or mechanism, as in 
those described above in this Section. The rotational dynamics 
of the nodal line arising from the interference of wave packets 
in the complex plane thus off'ers an elegant method to compute 
the lifetime for the interference features observed on the real 
axis. 

We consider two cases of head-on collision of two Gaus- 
sian wave packets, which depend on the relative magnitude be- 
tween the propagation velocity and the spreading rate of the 
wave packets | 30]. In the case where the relative propagation 
velocity is larger than the spreading rate of the wave packets, an 
average wrapping time is calculated to provide a lifetime for the 
interference features, while the rotational dynamics of the nodal 
line in the complex plane explains the transient appearance of 
the interference features observed on the real axis. In the case 
where the relative propagation velocity is approximately equal 
to or smaller than the spreading rate of the wave packets, the 
rapid spreading of the wave packets leads to a distortion of 
quantum caves. The infinite survival of interference features 
in this case implies that the wrapping time becomes infinity. 
However, the rotational dynamics of the nodal line clearly ex- 



plains the persistent interference features observed on the real 
axis. In both cases, the interference features are observed on 
the real axis when the nodal line is near the real axis; therefore, 
in contrast to conventional quantum mechanics, the rotational 
dynamics of the nodal line in the complex plane provides a fun- 
damental interpretation of quantum interference. 

The organization of this work is as follows. To be self- 
contained, in Sec. [2] we briefly describe the complex quantum 
Hamilton- Jacobi formalism as well as the QMF local structures 
and its associated PVF near characteristic points. In Sec.O first 
we present theoretical analysis of the Gaussian wave-packet 
head-on collision on the real axis and in the complex plane, 
and then demonstrate the interpretation of quantum interference 
with two cases in the framework of the complex quantum tra- 
jectory method. Finally, we present a summary and discussion 
in Seed 

2. Theoretical formulation 

2.1. Complex quantum Hamilton- Jacobi formalism 

Substituting the complex-valued wave function expressed 
by ^(x, t) = exp [iS (x, t)/h] into the time-dependent Schrodinger 
equation, we obtain the quantum Hamilton- Jacobi equation (QHJE) 
in the complex quantum Hamilton- Jacobi formalism. 



- — - — (— f -h V(x) -h — — 
dt 2m \dx I 2mi dx^ ' 



(1) 



where S{xj) is the complex action and the last term is the 
complex quantum potential. In real space, the QMF is de- 
fined by p{x, t) = dS (x, t)/dx, which, within the complex quan- 
tum Hamilton- Jacobi formalism, is analytically continued to 
the complex plane by extending the real variable x to a com- 
plex variable z = x -\- iy. The same is done with other relevant 
functions, such as the wave function and the potential energy. 
Quantum trajectories in complex space are then determined by 
solving the guidance equation 



dz 
dt 



p(z, t) 



(2) 



where time remains real- valued and the (complex) QMF is ex- 
pressed in terms of the (complex) wave function (through the 
complex action) as 



I T(z, oz 



(3) 



From this equation, nodes in the wave function correspond to 
QMF poles. Moreover, those points where the first derivative of 
the wave function vanishes will correspond to QMF stagnation 
points. 

2.2. Local structures of the quantum momentum function and 
its Polya vector field 

The dynamics of complex quantum trajectories is guided by 
the wave function through the QMF. The complex- valued QMF 
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Table 1 : Local structures or streamlines for the QMF and the associated PVF 
near a stagnation point or a pole. 



Taylor series around a stagnation point 





Stagnation point 


Pole 


QMF 


Spirals, circles, 


East- West and North-South 


or straight lines 


opening hyperbolic flow 


PVF 


Rectangular 
hyperbolic flow 


Circular flow 



can be regarded as a vector field in the complex plane. The tra- 
jectory dynamics is significantly influenced by both the QMF 
stagnation points and poles. Streamlines near a QMF stagnation 
point may spiral into or out of it, or they may become circles 
or straight lines. The QMF near a pole displays East- West and 
North-South hyperbolic structure 1671 [68|]. These local struc- 
tures around these characteristic points have also been observed 
for the one-dimensional stationary scattering problems includ- 
ing the Eckart and the hyperbolic tangent barriers ||47[|48|] . 

The PVF of a complex vector field f(z) is defined by its 
complex conjugate f*(z) I69l-[7lll. Thus, the PVF associated 
with the QMF, p(z, t) = px -\- ipy, is given by the vector field 
P(z, t) = px - ipy = (Px, ~Py)- This new vector field provides 
a simple geometrical and physical interpretation for complex 
circulation integrals 



p{z)dz = (^V-di + i(^V-dn, 



(4) 



where C denotes a simple closed curve in the complex plane, 
d£ = (dx, dy) is tangent to C and dn = (dy, -dx) is normal to C 
and pointing outwards. The real part of the integral in Eq. Q 
gives the total amount of work done in moving a particle along 
a closed contour C subjected to P, while its imaginary part gives 
the total flux of the vector field across the closed contour [ 69]. 

In Bohmian mechanics, quantum vortices form around nodes 
in two or more real coordinates ||26l 172148111 . Analogously, in 
the complex plane quantum vortices form around nodes of the 
wave function, the quantized circulation integral arising from 



the discontinuity in the real part of the complex action 11671 16811 . 
The PVF of a complex function contains exactly the same in- 
formation as the complex function itself, but it is introduced to 
interpret the circulation integral in terms of the work and flux 
of its PVF along the contour. Moreover, the PVF is the tangent 
vector field of contours for the complex-extended Born proba- 
bility density is^^lssll. Streamlines near a PVF stagnation point 
display rectangular hyperbolic structure, while streamlines near 
a PVF pole become circles enclosing the pole ||67l [68|1 . Local 
structures or streamlines for the QMF and its associated PVF 
are summarized in Table [T] 



2.3. Approximate quantum trajectories around a stagnation 
point 

Since the QMF is generally time-dependent, we can deter- 
mine approximate complex quantum trajectories around a stag- 
nation point (zo, ^o) in spacetime. We expand the QMF in a 
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where p{zo, fo) = has been used and the partial derivatives are 
evaluated at the stagnation point. Substituting this equation into 
Eq. ©, we obtain a first-order nonautonomous complex ordi- 
nary diff"erential equation for approximate quantum trajectories 



— = -{az+pt), 
at m 



(6) 



where (dp/dz)o = o; and (dp/dt)o = JS have been used and the 
origin has been moved to the stagnation point. General so- 
lutions for the nonautonomous linear diff'erential equation are 
given by H 



z(t) = e^'^^\z(0) + 



Jn Jo 



-as/m 



sds 



V \ ml 



(7) 



where z(0) is the starting point of a local approximate quan- 
tum trajectory. If we consider complex quantum streamlines at 
a specific time ts, the QMF p{z, ts) does not depend on time. 
Hence, the partial derivative of the QMF with respect to time 
is equal to zero, (dp/dt)o =13 = 0. Thus, the general solution 
in Eq. (O gives the approximate quantum streamlines near the 
stagnation point. 

2.4. Divergence and vorticity of the quantum momentum field 
and its Polya vector field 
The QMF first derivative contains the information about the 
divergence and vorticity of the quantum fluid in the complex 
plane |661. This can be shown as follows. The QMF diver- 
gence, which describes the local expansion/contraction of the 
quantum fluid, is given by 



OX oy 



(8) 



Analogously, the vorticity describing the local rotation of the 
quantum fluid is defined by the QMF curl. 



OX oy 



(9) 



Since the QMF is analytically extended to the complex plane, 
we use the Cauchy-Riemann equations to write the QMF first 
derivative as 



dp dpx ^ .dpy dpy dpx 1 ^ 
— = — — -h 1-— = — I—— = - (F -h iLl) . 

oz ox ox oy oy 2 



(10) 



Thus, the real and imaginary parts of the QMF first derivative 
determine the divergence and vorticity, respectively. Moreover, 
the complex quantum potential in E(L_ffl) can be expressed in 
terms of divergence and vorticity by ll66ll 
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(11) 



Additionally, we can also evaluate both the PVF divergence 
and vorticity, which are given by 

rp = V-P=^-^=0, (12) 
ox oy 

Qp = |VxP| = -^-^ =0, (13) 

OX oy 

respectively, where the Cauchy-Riemann equations for the QMF 
have been used. The vanishing divergence and vorticity indicate 
that the PVF associated with the QMF describes an incompress- 
ible and irrotational flow in the complex plane except at nodes 
of the wave function. Actually, this result follows from the fact 
that the PVF of a complex function /(z) is divergence-free and 
curl-free if /(z) is analytic, and vice versa |69]. 

The quantized circulation integral around quantum vortices 
originates from the work term of the PVF in Eq. (|4]) (671], 



r = Piz)dz = ^F'd£ = Innh. 



(14) 



The PVF near a pole is expressed in terms of polar coordinates 
by P = (y/27ir)ee, where y is the circulation and r is the radial 
distance from the center of the vortex. From Eq. ([13]), the PVF 
vorticity is zero everywhere except at poles. The PVF velocity 
near a pole varies inversely as the distance r from the core of the 
vortex. The circulation integral along a closed path enclosing 
the vortex is equal to y = Innh, which is independent of r. Al- 
though the quantum fluid described by the streamlines around 
a pole moves along a circular path, its vorticity is zero. These 
features indicate that the quantum vortex described by the PVF 
is a free or irrotational vortex [85,1 . 



2.5. Divergence and vorticity around a pole 

For a wave function with an ^-th order node at z ■ 



(z - ZpTfiz), we can evaluate the QMF first derivative, which 
yields 

dp ^ nih dps 

where Ps(z) is the smooth part of the QMF. Thus, the QMF first 
derivative can be approximated by the first term in the vicinity 
of a pole. For simplicity, we move the origin to the pole. Sepa- 
rating the first term in Eq. ([T5]) into its real and imaginary parts, 
we obtain the divergence and vorticity around the pole through 
Eq. dlOl), 



T = nh 



Q = nh 



4xy 



(jc2-h/)2' 

2(x^-/) 
(x2-F>;2)2' 



(16) 
(17) 



respectively. These approximate forms describe the local be- 
havior of the divergence and vorticity in the vicinity of the pole. 

Figure[T]presents the variations of both the QMF divergence 
and vorticity along the approximate streamlines in the vicinity 
of a pole. From now on, h = m = I and atomic units will be 
used throughout this work. In this case, we consider a wave 
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Figure 1 : (Color online) (a) The QMF displays hyperbolic flow around a pole. 
Particles start at the initial positions shown as dots and travel along streamlines 
from t = -n/3 to t = n/3. Variations of the QMF divergence and vorticity 
along: (b) streamlines 1 and 3 and (c) streamlines 2 and 4. 



function with one first-order node at the origin (n = 1). As 
discussed in Sec. 12.21 in Fig. [Ha) we observe how the QMF 
displays a hyperbolic flow around the pole. In addition, stream- 
lines 1 and 3 are parametrized by (x = +0. 1 sec t,y = ±0. 1 tan t) 
and streamlines 2 and 4 are parametrized by (x = ±0. 1 tan t,y = 
=f0.1 secO- In Figs. [Hb) andlUc), we show the variations of 
both the QMF divergence and vorticity given by Eqs. ([T6l) and 
([TT]) , respectively, along the streamlines shown in Fig.lTJa) from 
t = -71 / 3 to t = 7r/3. As can be noticed in these figures, the 
positive QMF divergence describes the local expansion of the 
quantum fluid when particles approach the pole. When parti- 
cles approach turning points near the pole, the QMF divergence 
vanishes; when particles leave the pole, the negative QMF di- 
vergence indicates the local contraction of the quantum fluid. 
On the other hand, the QMF vorticity describes the local ro- 
tation of the quantum fluid. During the whole process, it can 
be noticed that streamlines 1 and 3 display local counterclock- 
wise rotation, while streamlines 2 and 4 display local clockwise 
rotation. The QMF vorticity attains the extrema at the turning 
points. Therefore, when particles approach the pole, they re- 
bound as they experience a repulsive force. 
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3. Quantum interference 

3.1. Head-on collision of two Gaussian wave packets 
3.1.1. On the real axis 

We consider the Gaussian wave-packet head-on colHsion 
which, despite its simpHcity, can be considered as representa- 
tive of other more compHcated, reaUstic processes character- 
ized by interference (e.g., scattering problems, diffraction by 
sHts, etc.). This process can be described by the following total 
wave function 



^(x, t) = i//l(x, t) + i//r(x, t). 



(18) 



Each wave packet (L or R, left or right, respectively) is repre- 
sented by a free Gaussian function as 



<A(X, t) = At ^-(^-^^)V4^.^ro+/p(x-x, 



)/h+iEt/h 



(19) 



where, for each component. At = (2nd-j) and the complex 
time-dependent spreading is 



d-t = cro\l-\- 



iht 
2mcr^ 



(20) 



with the initial spreading ctq. From Eq. (l20l) , the spreading of 
this wave packet at time t is 



1 -h 



ht 

2mcr^ 



(21) 



Due to the free motion, Xt = xq -\- Vpt (Vp = p/mis the propa- 
gation velocity) and E = p^/2m, i.e., the centroid of the wave 
packet moves along a classical rectilinear trajectory. This does 
not mean, however, that the average value of the wave-packet 
energy is equal to E since, due to the quantum potential, the 
average energy is given hy E = p^/2m -h h^/Smo-Q [30]. We ob- 
serve two contributions in this expression. The former is associ- 
ated with the translation of the wave packet traveling at velocity 
Vp, while the latter is related to its spreading at velocity and 
has, therefore, a purely quantum-mechanical origin. The effec- 
tive momentum ps can then be written as ps = h/2crQ, which 
resembles Heisenberg's uncertainty relation. 

The relationship between Vp and plays an important role 
in effects that can be observed when dealing with wave packet 
superpositions |30]. Defining the timescale r = 2mcr^/h, we 
note that, if ^ «c r, the width of the wave packet remains basi- 
cally constant with time, cr^ ^ ctq (i.e., for practical purposes, 
it is roughly time-independent up to time t). This condition 
is equivalent to having an initial wave packet prepared with 

<^ Vp. Thus, the translational motion will be much faster 
than the spreading of the wave packet. On the contrary, if ^ » r, 
which is equivalent to » v^, the width of the wave packet in- 
creases linearly with time (cr^ ^ ht/2mcrQ), and the wave packet 
spreads very rapidly in comparison with its advance along x. Of 
course, in between, there is a smooth transition; from Eq. (|2T1) , 
it is shown that the progressive increase of cr^ describes a hyper- 
bola when this magnitude is plotted time. Thus, we can con- 
trol in the interference process the spreading and translational 
motions which determine the wave packet dynamics 



3.1.2. In the complex plane 

In conventional quantum mechanics, the interference pat- 
tern transiently observed on the real axis is attributed to the 
constructive and destructive interference between two counter- 
propagating components of the total wave function in Eq. ([TS]) . 
In contrast, in the framework of the complex quantum Hamilton- 
Jacobi formalism, the total wave function is analytically contin- 
ued to the complex plane from ^(x, t) to ^(z, 0- Therefore, two 
propagating wave packets always interfere with each other in 
the complex plane, and this leads to a persistent pattern (line) 
of nodes and stagnation points which rotates counterclockwise 
with time |60, 66]. 

The nodal positions in the complex plane can be determined 
analytically by solving the equation ^(z, t) = 0, resulting 



Zn(t) : 



i7i(n-\- 1/2) 



[imvp/h - (xo -Vpt) I (2croa^?)] 



(22) 



where n = 0, ±1, ±2, . . .. Here, we have assumed xr = -xl = 
xo and VpL = VpR = Vp. Splitting this expression into its real and 
imaginary parts, i.e., z«(0 = Xn(t) + iynO), we obtain 



Xn(t) =7r ^-h - - 

2 m 



Xot + VpT^- 



yn(t) = 7i\n-\- - 



2crl (vpt - xo) 



4 + 



(23) 
(24) 



respectively, where r = 2mcr^/h is the timescale for the Gaus- 
sian wave packet. In addition, dividing Eq. (l24l) by Eq. ([23l) 
yields the analytical expression for the (time-dependent) polar 
angle describing the angular position of the nodal line with re- 
spect to the positive real axis. 



6(t) = (tan)" 



yn(t) 



Xnit) 



= (tan) 



-1 



- {vpt - Xo) 



Xot + VpT^ 



(25) 



which does not depend on n. From this expression, we can 
calculate the rotation rate of the nodal line in the complex plane. 



oj(t) = 



dOit) 
dt 



h 



2mcr^ 



(26) 



where cr^ is given in Eq. ([2T]) . This equation indicates that the 
rotation rate is completely determined by the initial spreading 
of the Gaussian wave packet ctq. In addition, this rate is always 
positive and decays monotonically to zero as t goes to oo. From 
Eqs. ([23]) and (l24l) . we can also determine the node separation 
distance between two consecutive nodes 



d(t) = ^J[Xn+l(t) - Xn(t)]^ + [yn+l(t) " J«(0]^ 

nhcTt 



(27) 



2^2 



where ps = h/2crQ is the effective momentum. This distance is 
independent of n and it increases with time. Moreover, elim- 
inating t in Eqs. (l23l) and ([24b yields the nth node trajectory 
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describing the time evolution of the node given by 



yn = \^jxn-(2n^l)\—^\. 



ncr: 



(28) 



Consequently, these nodal trajectories with the same slope and 
different intercepts are parallel to each other. 

When t = 0, the initial angle of the nodal line is = 
(tan)"^(-xo/VpT). As described by Eq. (l26b , the positive ro- 
tational rate indicates that the polar angle of the nodal line in- 
creases monotonically with time. The nodal line thus rotates 
counterclockwise from the initial angle to a maximum or lim- 
iting value 600 = (tan)"^(VpT/xo) when t ^ 00. The angular 
displacement from ^ = to ^ = 00 is always equal to = 
Ooo - Oq = 71/ 2, because the product of the slopes of the nodal 
lines is equal to (tan ^0) (tan ^00) = -1. In particular, if both 
wave packets are initially very far apart (i.e., xq 00), but 
move with a finite velocity v, or they are at an arbitrary finite 
distance, but v = 0, the nodal line ends up aligned with the real 
axis. Otherwise, the nodal line starts at some angle ^0 and then 
evolves with the angular displacement AO = 71/ 2 until it reaches 
the limit angle ^00. In addition, the initial nodal line is perpen- 
dicular to all nodal trajectories in Eq. (1281) . Then, the nodal line 
rotates counterclockwise with time and it becomes parallel to 
nodal trajectories when t approaches infinity. 

Interference features are observed on the real axis only when 
the nodal line is near the real axis. When the nodal line co- 
incides with the real axis, the maximum interference features 
can be observed in real space. At this time, setting = in 
Eq. (|28]) , we recover the expression for the positions of nodes on 
the real axis, Xn = (n-\-l /2)A/2 where A = 27rh/mVp. During the 
evolution of the nodal line, its intersections with nodal trajecto- 
ries determine the positions of the nodes. As can be noticed, the 
time evolution of the nodal line in the complex plane therefore 
provides an elegant interpretation of quantum interference. 

3.2. Case 1: Vp > Vg 

We first consider the case where the relative propagation 
velocity is larger than the spreading rate of the wave packets. 
The following initial conditions for Gaussian wave packets are 
used: xql = -10 = -xqr, Vpi = 2 = -VpR and (Xo = ^^2. With 
these conditions, maximal interference occurs at ^ = 5 on the 
real axis and the propagation and spreading velocities are given 
by = 2 and = V2/4, respectively. 

3.2.1. Quantum caves with quantum trajectories 

Figure [2];a) displays complex quantum trajectories with the 
quantum caves consisting of the isosurfaces |^(z, 01 = 0.053 
and \d'^{zj)ldz\ = 0.106 from ^ = to ^ = 10 in a time- 
dependent three-dimensional Argand plot. As discussed in Sec. 12. 2 
the QMF local structures near stagnation points and poles pro- 
vide a qualitative description of the behavior of these trajec- 
tories. It is clearly seen that stagnation and vortical tubes al- 
ternate with each other, and the centers of the tubes are stag- 
nation and vortical curves. Trajectories display helical wrap- 
ping around the stagnation tubes and they are deflected by the 




Figure 2: (Color online) (a) Quantum caves and complex quantum trajecto- 
ries for the head-on collision of two Gaussian wave packets with the relative 
velocity larger than the spreading rate. These trajectories launched from the 
isochrone arrive at the real axis at t = 5. The caves are formed with the iso- 
surfaces mz,t)\ = 0.053 (pink/lighter gray surface) and \d'i'(z,t)/dz\ = 0.106 
(violet/darker gray surface), (b) Complex quantum trajectories displaying he- 
lical wrapping around the stagnation tubes (violet/darker gray) and hyperbolic 
deflection around the vortical tubes (pink/lighter gray). 



vortical tubes to show hyperbolic indentations in three dimen- 
sional space, as shown in Fig. Ob). These trajectories, which 
display complicated paths around stagnation and vortical tubes, 
depict how probability flows. In addition, trajectories from dif- 
ferent launch points wrap around the same stagnation curve and 
remain trapped for a certain time interval. Then, they sepa- 
rate from the stagnation curve. This counterclockwise circula- 
tion of trajectories can be viewed as a resonance process. This 
phenomenon characterizes long-range correlation among tra- 
jectories arising from difl'erent starting points. Trajectories may 
wrap around the same stagnation curve with difl'erent wrap- 
ping times and numbers of loops. Trajectories starting from 
the isochrone with the small initial separation may wrap around 
difl'erent stagnation curves and then end up with large separa- 
tions. In this way, interference leads to the formation of quan- 
tum caves and the topological structure displayed by complex 
quantum trajectories. 

3.2.2. Dislocations of the complex action 

The complex action function S (z, t) displays fascinating fea- 
tures in the complex plane. Decomposing this function into 
its real and imaginary parts S = S r -\- iS j, write the wave 
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Figure 3: (Color online) Real (a) and imaginary (b) part of the complex action 
for the complex-extended wave function associated to Eq. ( fTsl at t = 2.5. For 
the real part, the principal value of the multivalued phase function is shown. 
Blue to red degradation indicates the transition from lower to higher values of 
the corresponding functions (in real part of the complex action, the range goes 
from -n to n). 



function as i//(z, t) = Qxp(-Si/h) Qxp(iSR/h). According to this 
expression, the real and imaginary parts of the complex action 
determine the phase and amplitude of the wave function, re- 
spectively. Figure [3] displays the real and imaginary parts of 
the complex action for the complex-extended wave function in 
Eq. ([TSl) at ^ = 2.5. Figure Oa) displays the principal zone of 
the phase of the wave function in the range -n < diVg{S r) < n. 
Figure Ob) displays the imaginary part S j of the complex ac- 
tion in the complex plane. The vanishing of the wave function 
at nodes indicates that the imaginary part of the complex ac- 
tion tends to positive infinity at nodes. The peaks in Fig. Ob) 
correspond to the nodes in the wave function. 

The quantized circulation integral around a node in the wave 
function can be related to the change in the phase of the wave 
function and this integral can be expressed in terms of the PVF 
by 167J 

^ p(z)dz = (^V-di = 2nnh = AcSr, (29) 

where C denotes a simple closed curve in the complex plane. 
The quantized circulation integral around a node of the wave 
function originates from the discontinuity in its phase. The 
PVF displays counterclockwise circular flow in the vicinity of 



Figure 4: (Color online) (a) Complex quantum trajectory launched from z = 
-9.1 1016 - 1.17309/ and reaching the real axis at z = -0.3 when t = 5. As can 
be noticed, this trajectory undergoes a helical wrapping around an stagnation 
curve (blue/darker gray) and hyperbolic deflection around the vortical curves 
(red/light gray), (b) Vorticity and divergence along the trajectory. 



a node f67', ^68^ . Nodes in the complex-extended wave function 
in Eq. ([TS]) are first-order nodes (n = 1). As shown in Fig.Oa), 
if we travel around a first-order node counterclockwise along a 
closed path, it follows from Eq. (|29l ) that the phase displays a 
sharp discontinuity of 27i. Actually, the phase of the wave func- 
tion is a multivalued function in the complex plane. Thus, if 
we travel counterclockwise around a node, we then go through 
the branch cut from one Riemann sheet to another. Through a 
continuous closed circuit around a node continuing on all the 
sheet, the phase function generates a helicoid along the vertical 
axis. Analogous to the case in Bohmian mechanics, phase sin- 
gularities at nodes in the complex plane can be interpreted as 
wave dislocations 10, H^. 

3. 2. 3. QMF divergence and vorticity 

FigurelUa) shows a trajectory launched from z = -9. 11016- 
1.17309/ which later arrives at z = -0.3, when maximal in- 
terference occurs at ^ = 5, and Fig. lUb) presents the diver- 
gence and vorticity along this trajectory. When the particle ap- 
proaches a turning point, its velocity undergoes a rapid change 
and the divergence and vorticity display sharp fluctuations. When 
the particle approaches the vortical curve at position a along 
the direction of streamline 2 shown in Fig. [Ha), the trajectory 
displays hyperbolic deflection and the divergence and vorticity 
display analogous variations, as shown in Fig. [He). Around 
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Figure 5: (Color online) (a) QMF first derivative at z = 0. (b) The exact trajec- 
tories (thick curves) starting from the isochrone arrive at the real axis at f = 5 
with the approximate trajectories (thin curves) determined by Eq. around 
the stagnation curve along z = 0. 



position a, when the particle approaches the vortical curve, the 
positive divergence indicates the local expansion of the quan- 
tum fluid. When the particle arrives at the turning points, the 
divergence vanishes. Then, when it leaves the vortical curve, 
the negative divergence indicates the local contraction of the 
quantum fluid. On the other hand, the negative vorticity de- 
scribes the local clockwise rotation of the quantum fluid. 

When the trajectory displays helical wrapping around the 
stagnation curve, the particle is trapped between two vortical 
curves. When the particle approaches the turning points b and 
d along the direction of streamline 1 shown in Fig. [Ha) and ap- 
proaches the turning point c along the direction of streamline 
3, the divergence and vorticity in Fig. Ub) display analogous 
variations shown in Fig. [Hb). When the particle approaches 
and leaves the vortical curve, the divergence describes the local 
expansion and contraction of the quantum fluid in the vicinity 
of the vortical curve. The positive vorticity indicates the coun- 
terclockwise rotation of the quantum fluid. Finally, when the 
particle leaves the stagnation curve, it approaches the vortical 
curve at position e along the direction of streamline 4 shown in 
Fig-IHa). Again, the divergence and vorticity around position e 
display similar fluctuations shown in Fig.[TJc). Therefore, the 
variations of the divergence and vorticity around a pole ana- 
lyzed in Sec. 12.51 provide a qualitative description of the local 
behavior of the complex quantum trajectories in the vicinity of 
vortical tubes. 



3.2.4. Approximate quantum trajectories 

As discussed in Sec. 12.21 the QMF first derivative evaluated 
at stagnation points determines the local structure around these 
points 1 68]. Since it is noted from Eq. ([TSl) that there is always 
constructive interference at the origin, the origin is a stagnation 
point at all times. As an example. Fig. Oa) presents the QMF 
first derivative at the origin from t = Oiot = 10. The real part of 
the QMF first derivative indicates that the QMF displays con- 
vergent flow around the stagnation point until ^ = 4.851. Then, 
the real part of the QMF first derivative becomes positive and 
the QMF displays divergent flow. On the other hand, the QMF 
initially displays clockwise flow around the stagnation point. 
After t = 0.395, the QMF displays counterclockwise flow. As 
displayed in Fig.[2l trajectories exhibit helical wrapping around 
the stagnation curve along z = from approximately ^ = 3.5 to 
t = 6.5. During this time interval, the positive imaginary part 
of the QMF first derivative in Fig. [3a) describes the counter- 
clockwise flow of the quantum fluid. However, the real part of 
the QMF first derivative changes sign at r = 4.851. Particles 
initially converge to the stagnation point and they are gradually 
repelled by the stagnation point after t = 4.851. Finally, these 
particles depart from the stagnation point. Therefore, the QMF 
first derivative at the stagnation point qualitatively explains the 
behavior of trajectories in the vicinity of the stagnation point. 

Additionally, Fig. Ob) shows that the trajectories starting 
from the isochrone arrive at the real axis at ^ = 5 with the 
approximate trajectories around the stagnation point given in 
Eq. d?]). Here, the stagnation point in spacetime for Eq. d?]) is 
chosen to be z = and ^ = 5. The positions for the approximate 
trajectories are set to be the same as those for the exact trajecto- 
ries at ^ = 5. As shown in this figure, the trajectory determined 
by Eq. d?]) is a good approximation of the exact trajectory, pro- 
vided that the approximate trajectory is close to the stagnation 
point in spacetime. 

3.2.5. Wrapping time 

These complicated features for trajectories arise from the 
complex quantum potential in Eq. ([TT1) . Both the QMF diver- 
gence and vorticity characterize the turbulent flow of complex 
quantum trajectories. Moreover, the variation of the QMF vor- 
ticity provides a feasible method to define the wrapping time for 
a specific trajectory around a stagnation curve. The wrapping 
time can be defined by the time interval between the first and 
last minimum of Q comprising a region with the positive vortic- 
ity. Within this time interval, the positive vorticity describes the 
counterclockwise twist of the trajectory displaying the interfer- 
ence dynamics, and the stagnation points and poles significantly 
aff'ect the motion of the particle. For example. Fig. (Hb) in- 
dicates that the particle undergoes counterclockwise wrapping 
around the stagnation curve from t = 3.676 to ^ = 6.954 and 
hence the wrapping time is % = 3.278. Diff'erent trajectories 
have diff'erent wrapping times. Therefore, the average wrapping 
time can be utilized to define the 'lifetime' for the interference 
process observed on the real axis. 

In Fig. [6] we present the wrapping times that correspond to 
the first and last minimum of the vorticity, which comprise a re- 
gion with positive vorticity for trajectories launched from the 
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Figure 6: (Color online) Wrapping times corresponding to the first and last 
minimum of the vorticity, which comprise a region with positive vorticity for 
trajectories launched from the isochrone. These trajectories arrive at the real 
axis at f = 5, covering a range which goes from z(5) = -3.9 to z(5) = 3.9 with 
an increment of Az = 0.05. The vertical dotted lines denote the positions of the 
nodes on the real axis at f = 5. 



isochrone arriving at the real axis at ^ = 5. Since the pro- 
cess described in Eq. ([TSb in the complex plane is symmetric 
with respect to the origin, the data shown in Fig. [6l display 
the same symmetry. As displayed in Fig.[2tb), trajectories can 
wrap around stagnation curves with different wrapping times 
and numbers of loops. Figures [6] and [2tb) indicate that trajec- 
tories wrapping around stagnation curves with small rotational 
radius are trapped between vortical curves for a long time, and 
this leads to a long wrapping time for these trajectories. On the 
contrary, if trajectories wrapping around stagnation curves with 
large rotational radius, when they approach vortical curves, they 
experience a significant repelling force provided by QMF poles. 
Thus, these trajectories can easily escape capture by stagnation 
curves, and this results in a short wrapping time. In addition, it 
is noted in Fig. [6l that all trajectories display helical wrapping 
from approximately ^ = 4.3to^ = 5.9. Furthermore, the aver- 
age wrapping time for these trajectories is % = 3.24, and it can 
be used to define the lifetime for the interference process in this 
case. 

3.2.6. Rotational dynamics of the nodal line 

As described in Sec. 13. 11 two counter-propagating Gaussian 
wave packets interfere with each other at all times in the com- 
plex plane, and the nodal line rotates counterclockwise with 
respect to the origin as time progresses. Figure [3a) shows 
the evolution of stagnation points and nodes and nodal trajec- 
tories in the complex plane, and Fig. Ob) displays the time- 
dependent probability densities along the real axis. Initially, 
the interference of tails of two wave packets contributes to the 
string of stagnation points and nodes, and the initial nodal line 
with the angle = -51.34° is perpendicular to nodal trajecto- 
ries in Eq. ([28l) . Then, the nodal line rotates counterclockwise 
and lines up with the real axis at ^ = 5. At this time, the to- 
tal wave function displays maximal interference and the exact 
nodes form on the real axis. After ^ = 5, these two wave packets 
start to separate but keep interfering with each other in the com- 
plex plane, and the nodal line continues to rotate counterclock- 
wise away from the real axis. Finally, the angle of the nodal 
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Figure 7: (Color online) (a) Evolution of stagnation points and nodes of the 
wave function with the nodal line (black solid line) and the nodal line at t = oo 
(black dashed line). The arrows indicate the rotational direction of the nodal 
line. Nodal trajectories are shown as dotted lines passing through the nodal 
points. Thick black solid lines correspond to the nodal lines with 0(3.52) = 
-10° and 6(132) = +10°. Note that 600-60= ^/2. (b) Time-dependence of 
probability densities along the real axis. 



line tends to the limit angle Ooo = 38.66° when t approaches 
infinity, and the nodal line becomes parallel to nodal trajecto- 
ries. The rotational rate of the nodal line in Eq. ([26b decays 
monotonically to zero when t tends to infinity. In Fig.[7ta), the 
intersections of the nodal line and nodal trajectories determine 
the nodal positions, and the distance between nodes in Eq. dZTl) 
increases with time. Therefore, the interference process is de- 
scribed by the rotational dynamics of the nodal line with the 
angular displacement AO = 71/ 2. 

In conventional quantum mechanics, interference extrema 
transiently forming on the real axis are attributed to construc- 
tive and destructive interference between components of the to- 
tal wave function, as shown in Fig. [Ttb). In contrast, in the 
complex quantum Hamilton- Jacobi formalism, the interference 
features observed on the real axis are described by the counter- 
clockwise rotation rate of the nodal line in the complex plane. 
Since interference features are observed on the real axis only 
when the nodal line is near the real axis, we can define the life- 
time for the interference process corresponding to the time in- 
terval for the nodal line to rotate from 6 = -10° to ^ = -1-10°. In 
Fig. [3a), the lifetime for the interference features is A^ = 3.8. 
Therefore, compared to the conventional quantum mechanics, 
the complex quantum Hamilton- Jacobi formalism provides an 
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Figure 8: (Color online) Poly a trajectories displaying helical wrapping around 
the vortical tubes (pink/lighter gray surface) and hyperbolic deflection around 
the stagnation tubes (violet/darker gray sheets). 



elegant method to define the Hfetime for the interference fea- 
tures observed on the real axis. 

3.2. 7. Quantum caves with Poly a trajectories 

Although the QMF displays hyperbolic flow around a node, 
its associated PVF displays circular flow LSL .68] • Figure [8] 
shows that Poly a trajectories launched from the isochrone ar- 
rive at the real axis with quantum caves. In contrast to Fig.Oa), 
these trajectories display helical wrapping around the vortical 
tubes and hyperbolic deflection around the stagnation tubes. 
Both the PVF divergence and vorticity vanish everywhere ex- 
cept at poles; thus, trajectories display helical wrapping around 
irrotational vortical curves described by the PVF. 

3.3. Case 2: Vp < Vs 

Next, we consider the case where the relative propagation 
velocity is approximately equal to or smaller than the spreading 
rate of the wave packets. We use the following initial conditions 
for Gaussian wave packets: xql = -5 = -xqr, Vpi = I = -VpR 
and (To = V2/4. Maximal interference also occurs at ^ = 5 on 
the real axis and the propagation and spreading velocities are 
given by = 1 and = V2, respectively. 

3.3.1. Quantum caves with quantum trajectories and Poly a tra- 
jectories 

Figure [9] shows that quantum trajectories and Poly a trajec- 
tories starting from the isochrone reach the real axis at ^ = 5 
with quantum caves consisting of the isosurfaces of the wave 
function and its first derivative. Similar to the case shown in 
Fig. O quantum caves form around stagnation curves and vor- 
tical curves appearing alternately, but they are significantly dis- 
torted due to the rapid spreading of the wave packets. In addi- 
tion, quantum trajectories again display helical wrapping around 
the stagnation tubes and hyperbolic deflection near the vortical 
tubes. On the contrary. Poly a trajectories display hyperbolic de- 
flection near the stagnation tubes and helical wrapping around 
the vortical tubes. Again, trajectories launched from diff'erent 




Figure 9: (Color online) Quantum caves (a) with complex quantum trajectories 
and (b) with Polya trajectories for the head-on collision of two Gaussian wave 
packets with the relative velocity smaller than the spreading rate. These trajec- 
tories launched from the isochrone arrive at the real axis at r = 5. The caves 
are formed with the isosurfaces |^(z, 01 = 0.16 (pink/lighter gray surface) and 
\d^{z,t)ldz\ = 0.23 (violet/darker gray surface). 



Starting points show long-range correlation, and interference 
leads to the formation of quantum caves and produces compli- 
cated behavior of these trajectories. 

3.3.2. Rotational dynamics of the nodal line 

As shown in Fig. [2l quantum trajectories for Case 1 re- 
main trapped for a certain time interval between vortical curves, 
and then they depart from the stagnation curves. In contrast. 
Fig. [9ta) indicates that quantum trajectories can wrap around 
stagnation curves for an infinite time. This is in agreement with 
our previous statement that the wrapping time is a measure of 
the lifetime of the interference features; in this case, these fea- 
tures remain visible asymptotically and the wrapping time be- 
comes infinity. 

Figure [TOta) shows the evolution of stagnation points and 
nodes and nodal trajectories in the complex plane, and Fig.fTOtb) 
displays the time-dependent probability densities along the real 
axis. The nodal line starting with the initial angle ^0 = -87.14° 
rotates counterclockwise with respect to the origin and reaches 
the real axis at ^ = 5 when the maximal interference is observed 
on the real axis. Then, the nodal line rotates counterclockwise 
away from the real axis, and it approaches the limit nodal line 
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Figure 10: (Color online) (a) Evolution of stagnation points and nodes of the 
wave function with the nodal line (black solid line) and the nodal line at t = oo 
(black dashed line). The arrows indicate the rotational direction of the nodal 
line. Nodal trajectories are shown as dotted lines passing through the nodal 
points. Thick black solid lines correspond to the nodal lines with 0(1.09) = 
-10°. Note that 6oo - Oq = n/2. (b) Time-dependence of probability densities 
along the real axis. 



with the angle 6oo = 2.86° as t tends to infinity. Therefore, the 
nodal line rotates with the angular displacement = ;7r/2 to 
reach the limit nodal line parallel to nodal trajectories. 

When the nodal line is near the real axis, the interference 
features are clearly displayed on the real axis. As in Case 1, we 
can define the starting time of the interference process as the 
time for the nodal line with the angle 6 = -10°. In Fig.fTOta), 
^(1.09) = -10° in this case. In addition, the limit nodal line 
with the angle ^co = 2.86° is extremely close to the real axis. 
Therefore, the interference process starts at t = 1.09 and re- 
mains until t tends to infinity. Figure [T3b) indicates that the to- 
tal wave function starts to display the interference feature when 
t = 1.09 and the interference feature persists at long times. 

Figure [TT] presents the angle and the rotational rate of the 
nodal line in Eqs. ([25]) and (l26b for Case 1 and Case 2. As 
indicated in Eq. ([26b , the rotational rates for these two cases 
both decay monotonically to zero when t tends to infinity. For 
Case 1 , the angle of the nodal line alters relatively slowly and 
the rotational rate gradually decreases to zero. In contrast, for 
Case 2, the rotational rate shows a rapid decrease within the 
initial time interval in Fig. [TlT b). This dramatic change in the 
rotational rate reflects the fast rotation of the nodal line from 
the initial position to the vicinity of the real axis in Fig. [TlT a). 
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Figure 1 1 : (Color online) Rotational angle and rotational rate for the nodal line 
for Case 1 (vp > Vs) and Case 2 (vp < v^): (a) Angle of the nodal line in 
Eq. f25\ . (b) Rotational rate of the nodal line in Eq. ( |26l . 



4. Final discussion and concluding remarks 

In this study, quantum interference was explored in detail 
within the complex quantum Hamilton- Jacobi formalism. We 
reviewed local structures of the QMF and its associated PVF 
around stagnation points and poles, and derived the first-order 
equation for approximate quantum trajectories around stagna- 
tion points. Analysis of both the QMF divergence and vorticity 
along streamlines around a pole was employed to explain the 
complicated behavior of complex quantum trajectories around 
quantum caves. In addition, both the PVF divergence and vor- 
ticity vanish except at poles; hence, the PVF describes an in- 
compressible and irrotational flow in the complex plane. In 
contrast, both the QMF divergence and vorticity characterize 
the turbulent flow in the complex plane. 

We analyzed quantum interference using the head-on colli- 
sion of two Gaussian wave packets as an example. Exact de- 
tailed analysis was presented for the rotational dynamics of the 
nodal line on the real axis and in the complex plane. Complex 
quantum trajectories display helical wrapping around stagna- 
tion tubes and hyperbolic deflection around vortical tubes. In 
contrast. Poly a trajectories display hyperbolic deflection around 
stagnation tubes and helical wrapping around vortical tubes. 
For the case where the relative propagation velocity is larger 
than the spreading rate of the wave packets, during the inter- 
ference process, trajectories keep circulating around stagnation 
tubes as a resonant process and then escape as time progresses. 
Phase singularities in the complex plane can be regarded as 
wave dislocations. Then, the wrapping time for an individ- 
ual trajectory was determined by both the QMF divergence and 
vorticity, and the average wrapping time was calculated as one 
of the definitions for the lifetime of interference. For the case 
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where the relative propagation velocity is approximately equal 
to or smaller than the spreading rate of the wave packets, the 
distortion of quantum caves originates from the rapid spread- 
ing of the wave packets. Due to the rapid spreading rate, in- 
terference features also develop very rapidly, remaining visible 
asymptotically in time. The wrapping time becomes infinity, 
and this implies the infinite survival of such interference fea- 
tures. However, since the interference features are observed 
on the real axis only when the nodal line is near the real axis, 
the rotational dynamics of the nodal line in the complex plane 
off'ers a unified description to clearly explain transient or persis- 
tent appearance of the interference features observed on the real 
axis. Therefore, these results show that the complex quantum 
trajectory method provides a novel and insightful interpretation 
of quantum interference. 

The average wrapping time determined by both the QMF 
divergence and vorticity can be used as one of the definitions 
for the lifetime of interference. On the contrary, the PVF di- 
vergence and vorticity cannot be used to define the lifetime of 
interference because they vanish except at poles. However, the 
PVF of a complex function, such as the QMF, contains exactly 
the same information as the complex function itself [69]. Thus, 
it is suflScient to define the lifetime of interference as the av- 
erage wrapping time determined by the QMF divergence and 
vorticity. 

The head-on collision of two Gaussian wave packets with 
equal amplitudes was used as a model system to explore quan- 
tum interference in complex space. This problem is the proto- 
type of quantum systems displaying interference eff'ects, and it 
also exhibits basic features of quantum interference. A straight- 
forward generalization is to consider the head-on collision of 
two Gaussian wave packets with diff'erent amplitudes. As shown 
in this paper, these two wave packets interfere with each other 
in the complex plane for all times. Because of the diff'erent 
amplitudes for these two counter-propagating wave packets, we 
cannot observe an infinite number of nodes on the real axis. 
However, nodes can occur on the real axis at specific times. In 
this case, the nodal line displays not only rotational motion but 
also translational motion in the complex plane. Analogously, 
the dynamics of the nodal line in the complex plane clearly ex- 
plains the interference features observed on the real axis. A 
detailed analysis will be reported in our future studies. 

The current study concentrates mainly on quantum inter- 
ference arising from the head-on collision of wave packets in- 
volving no external potential. We have presented a compre- 
hensive exact analytical study for this problem, and these re- 
sults demonstrates that the complex quantum trajectory method 
can provide new physical insights for analyzing, interpreting, 
and understanding quantum mechanical problems. There are 
various quantum eff'ects resulting from quantum interference. 
Therefore, in the future, quantum interference incorporating in- 
teraction with external potentials during physical processes can 
be examined through analytical and computational approaches. 
Multidimensional problems displaying quantum interference de- 
serve further investigation within the complex quantum Hamilton- 
Jacobi formalism. 
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